System Of Nonlinear Equations Research Articles

HYSTERESIS OF IMPLICIT EQUATIONS IN HYPOPLASTICITY FOR SOIL MATERIALS WITH GRANULAR HARDNESS DEGRADATION

AbstractWe study a hypoplastic model for soil and granular materials stemming from geomechanical engineering which further incorporates effects of degradation of the granular hardness, therefore allowing for the description of environmental weathering. The governing system is described by a nonlinear system of transcendental-differential equations for stress and strain rate, which is investigated with respect to its long-time dynamic. Under deviatoric stress control, two different solutions of the underlying, implicit differential equations are constructed analytically. The spherical components of stress and strain rate converge asymptotically to an attractor and lead to the sparsification of material states. Whereas under cyclic loading-unloading carried out in a numerical simulation, finite ratcheting of the deviatoric strain rate is observed in the form of a square spiral.

Analog quantum simulation of partial differential equations

Abstract Quantum simulators were originally proposed for simulating one partial differential equation in particular -- Schrodinger's equation. Can quantum simulators also efficiently simulate other partial differential equations? While most computational methods for partial differential equations -- both classical and quantum -- are digital (they must be discretised first), partial differential equations have continuous degrees of freedom. This suggests that an analog representation can be more natural. While digital quantum degrees of freedom are usually described by qubits, the analog or continuous quantum degrees of freedom can be captured by qumodes. Based on a method called Schrodingerisation, we show how to directly map D-dimensional linear partial differential equations onto a (D+1)-qumode quantum system where analog or continuous-variable Hamiltonian simulation on D+1 qumodes can be used. This very simple methodology does not require one to discretise partial differential equations first, and it is not only applicable to linear partial differential equations but also to some nonlinear partial differential equations and systems of nonlinear ordinary differential equations. We show some examples using this method, including the Liouville equation, heat equation, Fokker-Planck equation, Black-Scholes equations, wave equation and Maxwell's equations. We also devise new protocols for linear partial differential equations with random coefficients, important in uncertainty quantification, where it is clear how the analog or continuous-variable framework is most natural. This also raises the possibility that some partial differential equations may be simulated directly on analog quantum systems by using Hamiltonians natural for those quantum systems.

Analytical solution to boundary layer flow and convective heat transfer for low Prandtl number fluids under the magnetic field effect over a flat plate

AbstractThe present study aims to quantify the flow field, flow velocity, and heat transfer features over a horizontal flat plate under the influence of an applied magnetic field, with a particular emphasis on low Prandtl number fluids. Nonlinear partial differential expressions can be incorporated into the ordinary differential framework with the use of appropriate transformations. This research utilizes the variational iteration method (VIM) to approximate solutions for the system of nonlinear differential equations that define the problem. The objective is to demonstrate superior flexibility and broader application of the VIM in addressing heat transfer issues, compared to alternative approaches. The results obtained from the VIM are compared with numerical solutions, revealing a significant level of accuracy in the approximation. The numerical findings strongly suggest that the VIM is effective in providing precise numerical solutions for nonlinear differential equations. The analysis includes an examination of the flow field, velocity, and temperature distribution across various parameters. The study found that improving temperature patterns, velocity distribution, and flow dynamics were all positively impacted by increasing the Prandtl numbers. As a result, this leads to the thickness of the boundary layer to decrease and improves heat transfer at the moving surface. Thus, the convection process becomes more efficient. When the strength of the magnetic field is increased, the velocity of the fluid decreases. This observation aligns with expectations since the magnetic field hampers the natural flow of convection. Notably, the convection process can be precisely controlled by carefully applying magnetic force.

Analytical Properties and Numerical Preservation of an Age-Group Sir Model: Application to the Diffusion of Information

Analytical Properties and Numerical Preservation of an Age-Group Sir Model: Application to the Diffusion of Information

Analyzing epidemics and designing policies with a modified SIR model

This research paper aims to study the impact of COVID-19 on the economy, where we have extended the standard epidemiological model, the Susceptible–Infected–Recovered (SIR) model. In the economy’s dynamics, we have incorporated aggregate demand and aggregate supply and further set up a system of three-dimensional nonlinear differential equations, which is more likely the Labor–Infected–Quarantine (LIQ) model, an extension of the SIR model. Here, we have derived the steady-state condition and analyzed the model’s stability. In addition, we have done a comparative static analysis of parameters with some policy implications. We aim to achieve the twin objective, i.e., controlling the pandemic and reviving the economy under the set of the policy mix.

Existence Results and Finite-Time Stability of a Fractional (p,q)-Integro-Difference System

In this article, we mainly generalize the results in the literature for a fractional q-difference equation. Our study considers a more comprehensive type of fractional p,q-difference system of nonlinear equations. By the Banach contraction mapping principle, we obtain a unique solution. By Krasnoselskii’s fixed-point theorem, we prove the existence of solutions. In addition, finite stability has been established too. The main results in the literature have been proven to be a particular corollary of our work.

Modified Newton–CAPRESB method for solving a class of systems of nonlinear equations with complex symmetric Jacobian matrices

In this paper, we contemplate addressing nonlinear problems involving complex symmetric Jacobian matrices. Firstly, we establish a parameter-free method called modified Newton–CAPRESB (MN–CAPRESB) method by harnessing the modified Newton method as the outer iteration and the CAPRESB (Chebyshev accelerated preconditioned square block) method as the inner iteration. Secondly, the local and semilocal convergence theorems of MN–CAPRESB method are proved under some conditions. Eventually, the numerical experiments of two kinds of complex nonlinear equations are presented to validate the feasibility of MN–CAPRESB method compared to other existing iteration methods.

Modified differential evolution to solve systems of nonlinear equations

Modified differential evolution to solve systems of nonlinear equations

A kinetic theory approach to modeling prey–predator ecosystems with expertise levels: analysis, simulations and stability considerations

A kinetic theory approach to modeling prey–predator ecosystems with expertise levels: analysis, simulations and stability considerations

Bioconvective flow of thermally radiative Casson nanofluid with aerobic microorganisms over a stretching surface

A mathematical model is proposed to analyze the flow characteristics of a thermally radiative, chemically reacting Casson nanofluid incorporating aerobic microorganisms over a stretching sheet. Nonlinear dimensionless ordinary differential equations are formulated through suitable similarity transformations. The resulting coupled system of nonlinear ordinary differential equations is numerically solved using the finite difference method with the bvp4c routine in MATLAB. The verification of the numerical solution’s accuracy is conducted by comparing it with previously established results under specific conditions pertinent to the present investigation. Skin friction coefficients, local Nusselt number, and Sherwood number are computed for various parameters and their implications for engineering applications are discussed. The investigation also explores the impact of different dimensionless parameters on velocity, temperature, mass concentration, motile microorganism density, and oxygen concentration profiles for both Casson and Newtonian nanofluids, offering a comprehensive analysis. The study highlights that suspending non-Newtonian nanofluids with fast-moving aerobic microorganisms and nanoparticles results in enhanced heat transfer rates, thus presenting potential benefits for practical applications such as the manufacturing of biofilms.

The null condition for systems of nonlinear wave equations and its related topics

In this article, we consider the Cauchy problem for systems of nonlinear wave equations, whose nonlinear terms depend mainly on derivatives of the unknown, with small, smooth, and compactly supported initial data. If the nonlinear terms have the critical power, we need some structural conditions to obtain the small data global existence. The null condition is one such condition, but recently some weaker conditions were also found. We discuss the small data global existence and the asymptotic behavior of solutions under these conditions. The corresponding result for nonlinear Schrödinger equations will be also discussed.

A well-balanced second-order finite volume approximation for a coupled system of granular flow

A well-balanced second-order finite volume scheme is proposed and analyzed for a 2×2 system of non-linear partial differential equations which describes the dynamics of growing sandpiles created by a vertical source on a flat, bounded rectangular table in multiple dimensions. To derive a second-order scheme, we combine a MUSCL type spatial reconstruction with strong stability preserving Runge-Kutta time stepping method. The resulting scheme is ensured to be well-balanced through a modified limiting approach that allows the scheme to revert to a well-balanced first-order scheme if it fails to capture or preserve the discrete steady state solution during its time evolution. Further, the scheme maintains the second-order accuracy away from the steady state. The well-balance property of the scheme is proven analytically in one dimension and demonstrated numerically in two dimensions. Additionally, numerical experiments reveal that the second-order scheme reduces finite time oscillations and provides better resolutions of the physical properties of the state variables, as compared to the existing first-order schemes of the literature.

Opers on the projective line, Wronskian relations, and the Bethe Ansatz

It is well-known that the spectra of the Gaudin model may be described in terms of solutions of the Bethe Ansatz equations. A conceptual explanation for the appearance of the Bethe Ansatz equations is provided by appropriate G-opers: G-connections on the projective line with extra structure. In fact, solutions of the Bethe Ansatz equations are parameterized by an enhanced version of opers called Miura opers; here, the opers appearing have only regular singularities. Moreover, this geometric approach to the spectra of the Gaudin model provides a well-known example of the geometric Langlands correspondence. Feigin, Frenkel, Rybnikov, and Toledano Laredo have introduced an inhomogeneous version of the Gaudin model; this model incorporates an additional twist factor, which is an element of the Lie algebra of G. They exhibited the Bethe Ansatz equations for this model and gave an interpretation of the spectra in terms of opers with an irregular singularity. In this paper, we consider a new geometric approach to the study of the spectra of the inhomogeneous Gaudin model in terms of a further enhancement of opers called twisted Miura-Plücker opers. This approach involves a certain system of nonlinear differential equations called the qq-system, which were previously studied in [20] in the context of the Bethe Ansatz. We show that there is a close relationship between solutions of the inhomogeneous Bethe Ansatz equations and polynomial solutions of the qq-system and use this fact to construct a bijection between the set of solutions of the inhomogeneous Bethe Ansatz equations and the set of nondegenerate twisted Miura-Plücker opers. We further prove that as long as certain combinatorial conditions are satisfied, nondegenerate twisted Miura-Plücker opers are in fact Miura opers.

Application of generalized Haar wavelet technique on simultaneous delay differential equations

In this paper, the Haar wavelet technique is applied to a system of simultaneous linear and nonlinear proportional delay differential equations. By implementing this method, the delay differential equations are transformed into a series of linear and nonlinear algebraic equations, respectively. These equations are then solved using numerical code constructed in MATLAB. The obtained approximate solutions are compared with exact solutions, demonstrating the accuracy and efficacy of the Haar wavelet method. Using the fixed point theorem, the existence and uniqueness of the solutions is established. The reliability and efficiency of the proposed technique are illustrated through two numerical examples. Additionally, a comprehensive error analysis is presented and the convergence result is discussed, providing a thorough examination of the robustness of the method. This study not only shows that the Haar wavelet technique is a powerful tool for solving delay differential equations but also sets the groundwork for future research in this area.

Energy estimated frequencies of standing-wave solutions to nonlinear Klein-Gordon systems in higher dimensions

A system of nonlinear elliptic equations is considered, where the nonlinear term is the sum of a quadratic form and a Sobolev sub-critical term. An extra assumption is introduced on the nonlinear term, which is minimal among the ones which guarantee the existence of standing-waves obtained by estimating frequencies of minimizing sequences with the energy functional.

Fractional view analysis of coupled Whitham- Broer-Kaup and Jaulent-Miodek equations

In this study, I systematically investigate the fractional Whitham-Broer-Kaup (WBK) system and the fractional Coupled Jaulent-Miodek (CJM) equation under the Caputo fractional calculus. The nonlinear fractional differential equation systems are investigated via the Aboodh transform iteration method and the Aboodh residual power series method, thus offering a thorough analytical investigation. The dynamics of the fractional WBK system are accomplished using the Aboodh transform iteration method, and the Aboodh residual power series method is employed to study the behavior of the fractional CJM equation. Using established solutions, we completely analyze their dynamics employing both symbolic computations and numerical simulations. Consequently, novel solutions are identified, and the behavior of these fractional systems in the Caputo operator sense are clarified. The results obtained show good agreement of convergence of the numerical and analytical solutions, highlighting the effectiveness of the schemes adopted for unraveling the complex dynamics of fractional nonlinear systems.

Modeling gas flow in low-permeability formations: an efficient combination of mixed finite elements and high order time integration schemes

Numerical simulation of gas flow in low permeability formations is a challenging task due to the high nonlinearity induced by (i) the compressibility of the gas, (ii) the Klinkenberg slippage effect and (iii) the Langmuir adsorption of the gas on the pore surface. Because of these nonlinearities, modeling gas flow in low permeability formations requires a great deal of computational effort. In this work, we develop an efficient numerical model using advanced spatial and temporal discretization methods for a simultaneous solution of the coupled equations of gas flow, cubic Peng-Robinson equation of state, slippage effect and Langmuir adsorption. The spatial discretization is performed with the lumped hybrid formulation of the mixed finite element method which is well adapted for fluid flow in heterogeneous porous media. The time integration is performed with high-order methods via the method of lines (MOL) which allows large time steps and efficient solution of the nonlinear system of equations. Numerical experiments, performed for gas extraction in a heterogeneous domain, point out the high efficiency of the new model since it can be until 10 times more efficient than the classical first-order time discretization method. Results of global sensitivity analysis show that the intrinsic gas-phase permeability and the Klinkenberg factor are the most influential parameters controlling the pumping rate in the case of gas extraction in a homogeneous domain.

A generalization of B-differentiable equations method for elastoplastic contact problems with dual mortar discretization

In this study, a novel method based on the B-differentiable equations (BDEs) is proposed to solve the three-dimensional (3D) static and dynamic elastoplastic contact problems with friction. The contact conditions are formulated as the BDEs system, which can precisely satisfy Coulomb’s law of friction. Based on the mortar segment-to-segment contact model, the contact constraints are enforced by dual Lagrange multipliers for nonconforming meshes of the contact boundary. The non-smooth and nonlinear system equations composed of equilibrium and contact equations are simultaneously solved by the B-differentiable damping Newton method (BDNM), which has global convergence property. In each iteration, firstly, some contact displacements and contact forces, called substituted contact variables, are expressed by others according to the different contact states and contact conditions. The equilibrium equation is then solved, where the stiffness matrix is updated by substituting the contact variables. At last, the substituted contact variables will be restored. In the proposed method, the contact flexibility matrix is no longer needed, which saves the time for calculating the contact flexibility matrix due to changes in the stiffness matrix resulting from the update of stresses in the iteration step. By comparing with the results obtained by finite element commercial software ANSYS, the accuracy of the proposed algorithm is verified, and the property of convergence and efficiency are demonstrated for the elastoplastic contact problems by several numerical examples.

Determining the parameters of the functioning for a nonlinear ballistic system in a real external environment

Purpose. Development of an approximate nonlinear model for the solution of the external ballistics problem with determination of the nonlinear system parameters and development of a methodology for calculating the aerodynamic characteristics of the cargo, located on the external suspension of unmanned aerial vehicle (UAV) in order to increase the efficiency of their delivery to the specified landing target by means of asymptotic approach with given parameters of the studied system and external load. Methodology. The development of an effective model was carried out using analytical and numerical research algorithms based on a nonlinear system of differential equations in a general form with time-varying coefficients. In order to obtain a solution to the nonlinear problem of external ballistics in a two-dimensional formulation, the assumption of a significant influence of the projection of the velocity function on the ordinate axis in relation to the component on the abscissa axis is introduced. The problem is reduced to the solution of a related system of differential equations with variable coefficients along the corresponding coordinates using the asymptotic approach for a small parameter of the coefficient of frontal aerodynamic resistance. Applied mathematical analysis and modeling have been used for the problem formulation considering studied environmental parameters. Findings. Analytical dependences of the nonlinear problem of ballistics and application of finite-element analysis (FEA) with respect to the cargo motion from the UAV in the presence of the initial speed and wind load in the plane of motion are proposed. It is shown that the obtained analytical solution is correlated with the direct numerical calculation of the basic differential equation with respect to the ordinate axis. Originality. A mathematical nonlinear model of the dynamic process is proposed, assuming the prevailing influence of the speed function along the ordinate axis compared to the function along the abscissa axis. To obtain an approximate analytical solution of the basic nonlinear system of differential equations with variable coefficients the asymptotic perturbation method is applied. The dependence for the axial displacement function is presented considering actual time-flight parameter. Practical value. The obtained analytical dependencies for estimating the time and distance reaching the target with the initial speed of movement and the presence of wind load can be used in applied problems of mathematical physics and engineering calculations of functional dependencies and control of the cargo delivery process and target reaching from an UAV. The obtained analytical results and the solution algorithm can be integrated into applied problems of mathematical physics and engineering calculations, particularly the development of ballistic system control algorithms.

Convex hull property for elliptic and parabolic systems of PDE

We study the convex hull property for systems of partial differential equations. This is a generalization of the maximum principle for a single equation. We show that the convex hull property holds for a class of elliptic and parabolic systems of non-linear partial differential equations. In particular, this includes the case of the parabolic p-Laplace system. The coupling conditions for coefficients are demonstrated to be optimal by means of respective counterexamples.